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Abstract 

In the present paper, given an evolving mixture of probability densities, 
we define a candidate diffusion process whose marginal law follows the same 
evolution. We derive as a particular case a stochastic differential equation 
(SDE) admitting a unique strong solution and whose density evolves as a mix- 
ture of Gaussian densities. We present an interesting result on the comparison 
between the instantaneous and the terminal correlation between the obtained 
process and its squared diffusion coefficient. As an application to mathematical 
finance, we construct diffusion processes whose marginal densities are mixtures 
of lognormal densities. We explain how such processes can be used to model 
the market smile phenomenon. We show that the lognormal mixture dynamics 
is the one-dimensional diffusion version of a suitable uncertain volatility model, 
and suitably reinterpret the earlier correlation result. We explore numerically 
the relationship between the future smile structures of both the diffusion and 
the uncertain volatility versions. 

Keywords: Stochastic Differential Equations, Mixtures of Densities, Mixtures of 
Gaussians, Mixtures of Lognormals, Risk-Neutral Valuation, Option Pricing, Volatility- 
Underlying Correlation, Smile Modeling. 

AMS classification codes: 60H10, 60J60, 91B28, 91B70 



* Draft. This working paper is downloadable at 



http : //www . damianobrigo . it 



D. Brigo: Mixture diffusion SDEs and volatility- asset covariance 



2 



Contents 

1 Introduction: SDEs and mixtures 

2 Diffusions whose densities follow mixtures of normal distributions [5] 



The smile phenomenon in option pricing 

3.1 The smile problem and market implied distributions . 

3.2 Local and stochastic volatility models 

3.3 Local volatility lognormal mixture diffusion dynamics 

3.4 Uncertain volatility geometric Brownian motion . . . 

3.5 Evolution of the volatility smile in the two models . . 



8 

8 

n 
12 

15 
17 



4 Conclusions and further research 



19 



D. Brigo: Mixture diffusion SDEs and volatility- asset covariance 



3 



1 Introduction: SDEs and mixtures 

Let us consider the scalar stochastic differential equation (SDE) 

dX t = f t (X t )dt + a t (X t )dW t , (1) 

of diffusion type, with deterministic initial condition X = x , and where {W t , t > 0} 
is a standard Brownian motion. We assume that 

(Al) The stochastic differential equation (OQ) characterized by the coefficients /, a, 
and by the initial condition x admits a unique strong solution, whose support 
is assumed to be the interval (b, oo) at all time instants. The symbol b denotes 
either a real number (typically 0) or — oo. 

Under (Al), we can analyze the distribution of our SDE's solution at all time instants. 
In describing the evolution of the distribution of a diffusion process, the Fokker- 
Planck partial differential equation is a fundamental tool. We assume that 

(A2) The unique solution X t of (pQ) admits a density p t that is absolutely continuous 
with respect to the Lebesgue measure in (b, oo) and that satisfies the Fokker- 
Planck equation: 

|^ = -^(M) + ^to), *(•) = ^(.). 

The other main ingredients in the present paper are mixtures of densities. More 
specifically, we will consider a basic parametric family of densities, say 

v = {p(-,9),9ee}, 

with open in R d , being d a suitable integer, and where all densities in the family 
share a common support (b, oo). We are interested in considering a particular mixture 
of densities in this family. In other words, we fix a set of non-negative weights 
Ai, . . . , A m , A > 0, J2i = 1- We take the space of all possible mixtures of densities 
in T> with fixed weights A: 

M(v, A) := {\ lP (; e 1 ) + ... + x mP (; e m ), e 1 ,...,e m e e}. 

We will refer to this set of densities as to the A-mixture family for T>, or shortly 
as to the "mixture family" when A is clear from the context. We are interested in 
finding an SDE whose solution X t has a density pt that follows a prescribed evolution 
in a given mixture family. More precisely, we require the curve t i— > p t , in the space 
of all densities, to match a given curve t i— > Y^i=\ \p{'-> ^i(^)) i n a given A4(V, A). 

Problem 1.1. Let be given a mixture family Ai(T>, A) of densities with support (b, oo), 
and a drift ft(x) satisfying 

(A3) lim f t {y)p{y,9) = for all t > 0, 9 e 6. 
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Let £(/, xq) denote the set of all real-valued diffusion coefficients o~f such that the 
related SDE (QP satisfies assumptions (Al) and (A2), and such that 

(A4) lim a({yfp{y, 0) = 0, Hm — (a{(y) 2 p(y, 9)) = for all t > 0, 9 £ 0. 

Assume Xq) to be non-empty. Then, given the curve t i— > X^Li ^*P(*> ^i(O) 
m M.(T), A) (where t \— > 9i(t) are C 1 -curves in the parameter space Q), find a dif- 
fusion coefficient in £(/, Xq) whose related SDE has a solution with density p t = 

This problem is the analogous of the drift search problem described and solved in 
Brigo and Mercurio (1998) and Brigo (2000). The solution of this problem is given 
by the following. 

Proposition 1.2. (Solution of Problem II. IT ) Assumptions and notation of Prob- 
lem li.il in force. Consider the stochastic differential equation 



dY t = f t (Y t )dt + a{(Y t )dW t , Y = x , (2) 

'A 
dt 



{4{y)f 1 



i=i 



E™i \p(y,9i(t)) 

+ / ft(x)y2xa>(x,6i(t))d 

Jb i=i 

Ij o* £ £(/, Xq), then the SDE solves Problem \l.l\ in that 

in 



±x ^ u R m) du)dx 



i=l 



Proof. Write the Fokker-Planck equation for the candidate diffusion Y with the can- 
didate solution YlT=i ^iPi'i ^iCO) already inside, and then back-out the diffusion co- 
efficient via two subsequent integrations starting from the lower point b. □ 

We will use a particular case of this proposition. 

Corollary 1.3. Assumptions as in the above proposition. If the basic densities 
p(-, 9i(t)) evolving in A4(T>, A) are the marginal densities of a family of (instrumental) 
SDEs 

dX\ = f t (Xt)dt + ai(Xi)dW t , x , p{x, 9i(t)) := p X i(x), 

all satisfying assumptions (Al), (A2), (A3) (with f l 's replacing f) and such that 
a 1 £ S(/*,a;o), then the solution of Problem U.l\) takes the form 
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An interesting particular case occurs when f satisfies 

m 

ft(y) = J2 A t(y)fi(y)- ( 4 ) 

i=i 

In this case the second term in the right hand side of (T3|) vanishes, and we have that 
requiring the marginal of our final SDE to be a given X-mixture of the marginals of 
the instrumental processes results in a drift and squared diffusion coefficients that are 
(state-dependent) A- "mixtures" of the drifts and squared diffusion coefficients of the 
instrumental processes. We move from the combinators X (density) to the combinators 
A (coefficients). 

Proof. It suffices to use the previous result by noticing that the Fokker-Planck equa- 
tion for the instrumental processes reads: 

dpjxAjt)) = 0[/*(s)pQr,0i(t))] + l d 2 [al(x) 2 p(^O i (t))] 
dt dx 2 dx 2 

and substituting in (T5]) the right hand sides of such Fokker-Planck equations. □ 

The problem solved in the above proposition could have been formulated to track 
a generic density-evolution t i— > q t which does not necessarily occur in a mixture 
family. In such a case, given assumptions analogous to (A3) and ( A4) , the appropriate 
diffusion coefficient would be 



Qt(y) 



^ dujdx+ / f t (x)q t (x)dx 



which generalizes the solution in the above proposition. We preferred to present 
directly the parametric case. In the next sections, we shall consider first an interesting 
application of the above corollary to the fundamental cases of mixtures of normal 
and lognormal families, and second an application of the latter to the option pricing 
problem in mathematical finance, in particular as a possible means to model the so 
called "smile" phenomenon. 

2 Diffusions whose densities follow mixtures of nor- 
mal distributions 

Mixtures of normals are ubiquitous in statistics, representing the standard in many 
applications. Also, mixture of normals are often used in econometrics to model time 
series that have tails fatter than the Gaussian. See for example Alexander (2001) for 
a discussion on normal mixtures applied to financial data, also in comparison with 
different distributions. In general, normal mixtures represent in a sense the least 
departure from the Gaussian family allowing for skewness and kurtosis different from 
the Gaussian ones. However, to the best of our knowledge, no explicit attempt has 
been made so far to design continuous time diffusion models displaying this kind of 
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marginal distributions. The subject can be relevant at least in mathematical finance, 
as we will show when addressing the smile problem, but bears also an interest of 
its own in the study of the interaction between a diffusion process dynamics and 
particular families of distributions. 

Let us then start from the normal family, that we parameterize via its first two 
moments. In this case 

= {PN(m,v2), m, v G R}, 6 = -oo. 

We are given a curve in the A-mixtures of normals, 1 1— > Y^T=i ^iPAf(mi(t),Vi(t) 2 ), and we 
wish to find a diffusion process compatible with such a law. It is easy to see that if 
we consider the instrumental processes 

dXi = fii(t) dt + a\t)dW t , [ ^(s)ds = m^t), f (a i (s)fds = v i (tf 

Jo Jo 

with deterministic \i and cr's and null initial condition, by applying the above corollary 
in the particular case (jl]) we end up with the diffusion process whose drift and squared 
diffusion coefficient are given respectively by 



m m 

fM = E Kiv)Atl °t(y) 2 = EA^ytt) 2 , (5) 

i=l i=l 

Now assume the coefficients t *—>■ /i(£)'s and t i— > cr(t)'s to be at least C 1 (and 
hence bounded on all finite time intervals), with the cr's bounded away from zero, 
&i{t) > L > for all % and t. There are possible problems for a regular behaviour of 
the above / and when (t, y) (0,0). In order to avoid this, we may decide to 
modify our coefficients by imposing //«(£) = fl, crj(t) = cr for all i and t G [0, e), with e 
a given positive real number, typically small. We can then assume suitable transitory 
trajectories for the /i's and the cr's in say [e, 2e) that recover the correct integrals m 
and v 2 for times larger than 2e. The desired mixture is thus unmatched only in the 
(typically negligible) time interval [0, 2e]. Now A|(y) = A, for alH < e and all y. 

From the above assumptions, since by definition < A\(y) < 1 for all t and y, 
we have also that / and cr^ are bounded in all regions with bounded time, implying 
non-explosion for the related SDE in all intervals of the kind t < T for some T > 0. 
Furthermore, / and (cr-^) 2 are also seen to be C 1 in both t and y. This in turn implies 
that / is locally Lipschitz, whereas the fact that the cr(t)'s are bounded away from 
zero implies that so is in both t and y, thus ensuring that (cr-^) 2 G C 1 implies 
<j3 g C 1 . Therefore a? is locally Lipschitz, and now we have all the elements to apply 
Theorem 12.1 in Section V.12 of Rogers and Williams (1996) to conclude that our 
SDE admits a unique strong solution. We have thus produced a diffusion process 
compatible with a given mixture of normal densities, and can state the following 
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Theorem 2.1. (SDEs whose marginal law follows a given normal mixture) 

Consider an SDE 

dY t = f t (Y t )dt + a({Y t )dW u Y = 0, 

with drift and diffusion coefficient f and o~* as in where the m 's and v 's are 
C 2 time functions. Assume moreover that in an initial time interval t e [0, e) we 
have rrii(t) = fit and Vi(t) 2 = o 2 t. Then the considered SDE admits a unique strong 
solution and its solution has as marginal density the X-normal mixture 

rn 

VYt = ^iPM{m i {t),v i {t)' 2 ) 
t=l 

at time t. 

An interesting feature of the obtained process concerns the covariance between 
the process itself and the diffusion coefficient in its dynamics. This quantity is of 
interest, for example, in mathematical finance. Let us denote by "corr t " the correla- 
tion between two random variables, and by "cov t " the covariance, both conditional 
on the information available at time t, the time being omitted if t = 0. We have, for 
the above SDE, the following "instantaneous" correlation between the instantaneous 
change in the process and the instantaneous change in the related diffusion coefficient 
at a given instant: 



sJd{Y) t d{a'(Y)), sfdYTdT, sj d a'(Y t ) d a'(Y,) 



as is obvious from the fact that the diffusion coefficient is a deterministic function of 
the current value of Y . 

However, things are rather different for the terminal correlation. A straightfor- 
ward if lengthy computation is needed to show that 

cov^,^) 2 ) = E (Y t a{(Y t ) 2 )-E (Y t ) E (a{(Y t ) 2 ) 

m / m \ / m N 

i=l \i=l / \i=l 

Consider now the case where all the means in the normal mixture densities 
are equal: //j(-) = /i(-) for all i, and correspondingly m^-) = m(-). In this case 
ft(y) = n(f) for all y and, perhaps surprisingly, especially if compared to the perfect 
instantaneous correlation, the above formula gives 

corr^,^) 2 ) = cov(Y t ,a?(Y t ) 2 ) = 0. 

This is a case where the instantaneous correlation is 1, whereas the terminal 
correlation after a time t, no matter how small, is 0. Thus we have a stochastic process 
whose instantaneous changes are perfectly correlated with the instantaneous changes 
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of its squared diffusion coefficient, whereas at any time its value has correlation 
with the squared diffusion coefficient value. It seems then that the squared diffusion 
coefficient has a special shape that immediately "decorrelates" itself from the process 
even after an arbitrarily small time. It is not difficult to prove an analogous statement 
for the average squared diffusion coefficient: 



Such interesting results on terminal correlation versus instantaneous correlation 
will be further discussed in the financial applications. 

Finally, going back to the above theorem, we notice that by reasoning along the 
same lines, diffusions displaying mixtures of lognormals as marginal densities can be 
easily obtained. Indeed, if Y is the diffusion process from the above theorem, we can 
set St '■— exp(Y t ) and derive easily the SDE for S via Ito's formula: 



implies that S is distributed as a mixture of lognormals if Y is distributed as a mixture 
of normals. 

3 The smile phenomenon in option pricing 

Now we briefly review a stylized version of the smile problem in financial modeling and 
explain the possible use of diffusions whose marginals follow A-mixtures of lognormals. 

3.1 The smile problem and market implied distributions 

Let us consider a financial market with a "money market account" process B t , with 
positive deterministic instantaneous interest rate r(t) > 0, so that dB t = r(t)B t dt. 
Let us also consider a process S t modeling the evolution of some traded financial 
(risky) asset in our market, typically a stock. 

The resulting financial market might admit arbitrage opportunities. A sufficient 
condition which ensures arbitrage-free dynamics is the existence of an equivalent mar- 
tingale measure Q, sometimes termed risk-neutral measure. An equivalent martingale 
measure is a probability measure that is equivalent to the initial one and under which 
the process {St/B t : t > 0} is a martingale. Let us assume, in line with the basic 
Black and Scholes (1973) setup, that the risk-neutral dynamics of S is modeled by 




dS t = 



S t ft(HSt))dt + \S t a({\n{S t )fdt + S t a{(\n(S t ))dW t , S = 1. 



It is a straightforward exercise to verify that the equality 



Ps t (y) =PY t Qny)/y 



dS t = r(t)S t dt + u(t)S t dW t , S = s , te [0, T], 



(6) 
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where so is a positive deterministic initial condition, and v is a well-behaving deter- 
ministic function of time (instantaneous volatility). The above process is a geometric 
Brownian motion and the probability density ps t of S t , at any time t, is lognormal. 
Indeed, 

ln^ ~tf(R(p,t) - \V{t)\ V{t) 2 ), R(a,t) := f r(s)ds, V{tf := [ u{s) 2 ds. (7) 

^0 ~ Ja JO 

When a = 0, we write shortly R(t) for R(0, t). In the above equation §6§ we modeled 
directly the risky asset dynamics under the unique equivalent martingale measure, so 
that W is assumed to be a Brownian motion under that measure, and it is immediate 
to check that S t /B t is indeed a martingale. It is this dynamics that matters when 
pricing options, as opposed to the real world one, which is related instead to historical 
estimation, statistical analysis and similar matters. Indeed, by applying the results by 
Harrison and Pliska (1981), the unique no-arbitrage price for a given Tt- measurable 
contingent claim H T e L 2 (Q) is V t = B t E Q {H T / ' B T \T t } =: B t E?{H T /B T } where 
\T t : t > 0} denotes the filtration associated to the process S. The contingent claim 
is said to be a simple one when it is of the form = h(ST) for a suitable function 
h. 

One of the most common simple claims is a European call option written on the 
stock, with maturity T and strike K, which pays H = (St — K) + at time T. Its 
price is obtained by computing the expectation of the discounted payoff according 
to the lognormal distribution implied by ([7]), leading to the celebrated Black and 
Scholes (1973) call option formula, which we denote by "BSCall" and whose explicit 
expression we omit for brevity: 

E$i(S T - K) + /B(T)} = BSCall(S , K, T, R(T), V(T)). 

The quantity V(T)/yT is the (average) volatility of the option, and according to 
this formulation, does not depend on the strike K of the option. Indeed, in this 
formulation, volatility is a characteristic of the stock 5* underlying the contract, and 
has nothing to do with the nature of the contract itself. In particular, it has nothing 
to do with the strike K of the option. 

Now take two different strikes K\ and K2. Suppose that the market provides 
us with the prices of two related options on our stock with the same maturity T: 
MKTCall(S ,^i,T) and MKTCall(So, K 2 ,T). 

Life would be simple if the market followed Black and Scholes' formula in a con- 
sistent way. But is this the case? Does there exist a single volatility parameter V(T) 
such that both the following equations hold? 

MKTCall(S , , K x , T) = BSCall(S , , K x , T, R(T), V[T)), 

MKTCall(So, K 2 , T) = BSCaU(Sb, K 2 , T, R(T), V(T)). 

The answer is a resounding "no" . In general, market option prices do not behave 
like this. What one sees when looking at the market is that two different volatilities 
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V(T, Ki) and V(T, K 2 ) are required to match the observed market prices if one is to 
use Black and Scholes' formula: 

MKTCall^o, K x , T) = BSCall(S , K x , T, R(T), V(T, K x )), 

MKTCall(So, K 2 , T) = BSCall(S , K 2 , T, R(T), V(T, K 2 )). 
In other terms, each market option price requires its own Black and Scholes (im- 



plied) volatility V MKT (T, K)/yT depending on the option strike K. 

The market therefore uses Black and Scholes' formula simply as a metric to express 
option prices as volatilities. The curve K 1— > v MKT (T, K) / \/T is the so called volatility 
smile of the T-maturity option. If Black and Scholes' model were consistent along 
different strikes, this curve would be flat, since volatility should not depend on the 
strike K. Instead, this curve is commonly seen to exhibit "smiley" or "skewed" 
shapes. 

Clearly, only some strikes K = Ki and maturities T = Tj are quoted by the mar- 
ket, so that usually the remaining points have to be determined through interpolation 
or through an alternative model. Interpolation in K, for a fixed maturity T, can be 
easy but it does not give any insight as to the underlying stock dynamics compatible 
with such prices. 

Indeed, suppose that we have a few market option prices for expiries T = Tj and 
for a set of strikes K = K^. 

For each fixed T = Tj, by smooth interpolation we can obtain the price for every 
other possible K, i.e. we can build a function K 1— > MKTCall(So, K, T). Now, if this 
strike- K price corresponds really to an expectation, we have 

poo 

MKTCall(S ,iT,T) = e" J ° r{s)ds \S T - K) + = e~'° < s)ds / (x - K)p T (x)dx, (8) 

where pt is the true risk-neutral density of the underlying stock at time T. If Black 
and Scholes' formula were consistent, this density would be the lognormal density, 
coming for example from a dynamics such as ([6]), i.e. px = Ps T - We have seen that 
this is not the case in the market. However, by differentiating the above integral 
twice with respect to K we see that, see also Breeden and Litzenberger (1978), 

so that by differentiating the interpolated-prices curve we can find the density pr 
of the underlying stock at time T that is compatible with the given interpolated 
prices. Nevertheless, the method of interpolation may interfere with the recovery of 
the density, since a second derivative of the interpolated curve is involved. Moreover, 
what kind of dynamics, alternative to (jUJ), do the densities Pt 1 iPt 2 i ■ ■ ■ ^Pt p ■ ■ ■ come 
from? 
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3.2 Local and stochastic volatility models 

A partial answer to these issues can be given the other way around, by starting from 
an alternative dynamics. Indeed, assume that 

dS t = r(t)S t dt + a(t,S t ) S t dW u S = s , (9) 

where a can be either a deterministic or a stochastic function of St- In the latter 
case we would be using a so called "stochastic- volatility model" , where for example 
<r(t, S) = with £ following a second stochastic differential equation, such as: 

dm 2 ) = Kt, m 2 )dt + x (t, m 2 )dz t , 

with the important specification 

dZ t dW t = p dt. 

Instead, in the so-called "local volatility models" the diffusion coefficient a(t, S t ) is a 
deterministic function of S t . 

One feature of stochastic volatility models that is usually deemed to render them 
superior with respect to local volatility models is "instantaneous decorrelation" . In- 
deed, for stochastic volatility models we can have 

Coii{dS t , da 2 {t,S t )) =p< 1 

whereas 

CoTv(dS t , da 2 (t,S t )) = 1 

for local volatility models, including our mixture dynamics models. This superiority 
no longer holds when considering terminal correlations, as we will remark later on. 
For the time being we concentrate on deterministic a(t, -)'s, leading to local- volatility 
models, such as for example cr(t,S) = r\ S 1 (CEV model, see Cox (1975)), where 7 
ranges in a suitable interval and where 77 is a positive deterministic constant. Below 
we will propose a new cr(i, •) of our own, flexible enough for practical purposes. 

We have seen above how the "true" risk- neutral densities px t , Pt 2 , ■ ■ ■ > PTj , ■ ■ ■ of 
the underlying asset are linked to market option prices through second-order differ- 
entiation. The problem we will face is finding a dynamics alternative to ([6]) and as 
compatible as possible with the densities pr 1 , Pt 2 > • • • 1 PTj , ■ ■ ■ ideally associated with 
market prices. This will be done by fitting directly the prices implied by our alterna- 
tive model to the market prices MKTCall(So, K, T) for the considered set of strikes 
Ki and maturities Tj. To further clarify this point, it may be helpful to explain 
explicitly how an alternative dynamics such as (Q leads to a volatility smile to be 
fitted to the market smile. The way in which an alternative local-volatility model 
dynamics generates a smile is clarified by the following stylized operational scheme: 

1. Set the pair (T, K) to a starting value; 
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2. Compute the model option price 

n(T, K) = e-fi r{s)ds E®(S T -K) + 
with S obtained through the no-arbitrage alternative dynamics 0. 

3. Invert Black and Scholes' formula for this strike and maturity, i.e. solve 

IT(T, K) = BSCall(S , K, T, R(T), V(T, K)) 
in V(T, K), thus obtaining the model implied volatility V(T, K). 

4. Change (T,K) and restart from point 2 until the last maturity/strike pair 
(T, K) is reached. 

The fact that the alternative dynamics is not lognormal implies that we obtain curves 
in the strike K > V(T, K) that are not flat. Clearly, one needs to choose a(t, •) 
flexible enough for the surface (T, K) i— > V(T, K) / yT to be able to resemble or even 
match the corresponding volatility surfaces coming from the market. Indeed, the 
model implied volatilities V(Tj,Ki)/ \/Tj corresponding to the observed strikes and 
maturities have to be made as close as possible to the corresponding market implied 
volatilities V MKT (Tj, Ki)/WTj, by acting on the coefficient u(-,S) in the alternative 
dynamics. 

At this point it should be clear why a A-mixture dynamics can be of help. Existing 
local volatility models have either too little flexibility to calibrate a large number of 
points in a volatility surface, or are specified in a too general way requiring interpo- 
lation and other possibly dangerous artifices in order to be implemented. The CEV 
model for example has only one more parameter 7 with respect to the basic Black 
and Scholes "flat" model in the time-homogeneous case, so that its fitting capabilities 
are rather poor. Dupire's (1997) approach is quite general, but if applied straight- 
forwardly in its most general form, it requires a continuum of traded strikes and 
maturities, with the possible interpolation problems observed above, not to mention 
the lack of guarantees on existence of solutions for the resulting SDE. 

3.3 Local volatility lognormal mixture diffusion dynamics 

Consider instead our approach and write the A-lognormal mixture diffusion. Set the 
instrumental processes to Black and Scholes processes, 

dX\ = r(t)X l t dt + Ui(t)XidW t , s 

and derive the diffusion coefficient corresponding to f(t,y) = r(t)y. Notice that 
here all instrumental processes have the same drift fa — f as the final process. In 
this particular case, taking into account the lognormal marginal distributions of the 
instrumental processes, by applying (j3J) under we have 

m 
i=l 
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Ht>y) 



Aj PAr(lnso+fl(t)-y»(i) 2 /2, K(t) 2 )(l n 2/) 



YJj=i \ PN{\ns +R(t)-v 3 (ty/2, v J (t)2)(^ny) , 



where R and V's are defined as in (J7J). In Brigo and Mercurio (2001b) we show that 
the SDE resulting from such coefficients, i.e. 



admits a unique strong solution, provided one takes suitable regularity conditions on 
the time functions t i— > z^(t)'s analogous to those on the t > cr l (t)'s illustrated in 
Section [2j Indeed, it is straightforward to prove such existence and uniqueness result 
starting from the proof given in Section [2] for the normal-mixture case. 

For our process S in ( ITOl . we confirm the curious result obtained in Section [2] 
from the comparison between instantaneous and terminal correlations. Notice that 
also in this case we have 



considered to be a drawback of local volatility models. Yet, when considering terminal 
correlations, things change considerably. 

Theorem 3.1. (Terminal correlation between underlying asset and average 
percentage variance in the lognormal mixture dynamics model for the 
smile) 

Consider the random variable 



dS t = r(t)S t dt + a mix (t, S t ) S t dW t , s 



(10) 



coiT t (dS t , da^Jt, S t )) = 1, 




v(T)/T being the "average percentage variance" of the process S. Then 



corr^a^^T, St), St) = 0, and corr(v(T), St) = for all T 



(11) 



Proof. First we show that 



corr(aL(T,S T ),S T )=0, 



i.e. that 




(12) 



e{oUt,St) s t } = / J2 A ^y)^ T ) 2 ypsr(y)dy = 




1=1 




^Ki{T,y)vi{Tfy ^ Xjp^^dy = / ^ A^(?/)z/ 4 (T) 2 ?/ dy 



m 



= s efo^J2^MT) 2 



i=i 
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given the definition of the A's. Similarly, one computes 

m 

E{alST,S T )} = J2*MT? 

i=i 

from which ffl~2l) follows. 



To show the other equality, notice that 

dv(t) = o^(t,S t )dt, 

and compute 

d(v(t)S t ) = S t crl ix (t, S t )dt + r(t)v(t)S t dt + (. . .)dW t . 

Taking expectations and Fubini's theorem 

dE(v(t)S t ) = E(S t al^t, S t ))dt + rE(v(t)S t )dt. 

Set At = E(S t (rl iix (t, S t )), which we computed above, and C t = E(v(t)S t ). The above 
equation reads 

C t = r(t)C t + A t , 

whose solution is 

Q f = e f r(s)ds / e -S;r{s)ds Adu _ 



By carrying out the computations one obtains E(v(t)S t ). At this point it is easy to 
prove that 

E(v(t)S t ) - E(v(t))E(S t ) = 
by computing E(v(t)) trough Fubini's theorem. □ 

The above result is partly weakened by the fact that correlation is not a satisfac- 
tory measure of dependence outside the Gaussian world. However, a striking feature 
remains of two processes that are instantaneously perfectly correlated but such that 
for any infinitesimal time T = e have zero terminal correlation. 

Let us now set apart this correlation result and go back to understanding the 
reason why a A-mixture dynamics is particularly appealing when pricing options. 
One of the main reasons lies in the price becoming a linear combination of prices 
with respect to underlying assets modeled according to the instrumental processes. 
Indeed, one has immediately, for a call option, 

U{T,K)=e-£^ ds E ( l{{S T -K)+}=e-ti^ ds / (y - K) + T >*Px' T (v)dy 

J o i=l 

m „oo m 

= J2ke~^ T r(s)ds / {y ^ K y px ^ y)dy = J2\ l BSCa\\(S ,K,T,R(T) } V l (T)) , 
i=i J o i=l 

the last equality following from the geometric Brownian motion structure of the 
underlying instrumental processes. This procedure is very general, and the price of 
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a European-style simple claim is always the linear combination of the corresponding 
prices for the instrumental processes. In the lognormal- mixture case, when pricing a 
call option we obtain a linear combination of Black and Scholes prices. This is very 
appreciated by traders, who usually prefer to contain departures from the lognormal 
distribution and the corresponding Black and Scholes formula. In a sense, when in 
need of generalizing a lognormal distribution, a mixture of lognormals is the least 
departure from the original lognormal paradigm. 

Important calibration benefits that should not go unnoticed are a consequence of 
the fact that a A mixture-of-lognormals dynamics can price call options analytically. 
This is very helpful for calibrating the model to the market. In such a case one runs 
an optimization to find the values of the parameters Vi and \ that best reproduce 
a given set of market prices, and the target function of this optimization can be 
computed in closed form without resorting to numerical methods such as Monte Carlo 
simulation, trees, and finite difference schemes. Notice also that since we are free to 
select an arbitrary number m of instrumental processes, in principle our diffusion 
model features a limitless number of calibrating parameters. Once the model has 
been calibrated, one can use it to price more complicated (for example early exercise 
or path dependent) claims that have no quoted price in the market. Monte Carlo 
simulation through the Euler or Milstein discretization schemes (see for example 
Kloden and Platen (1995)) applied to our dynamics or recombining trees in the spirit 
of Nelson and Ramaswamy (1990) can be attempted, thanks to the explicit diffusion 
dynamics we provided. 

In the present paper we have presented the mixture diffusion model in its most 
mathematical aspects. Other advantages and characteristics of these models and of 
their variants, based on shifted dynamics, and numerical investigation and calibra- 
tions to market data have been illustrated in Brigo and Mercurio (2000a, 2000b, 
2001a, 2001b). 

The introduction of general drift rates fii(t) not necessarily all equal to r in the 
instrumental processes and possible mixtures of densities coming from hyperbolic-sine 
processes are considered in Brigo, Mercurio and Sartorelli (2003). 

A study of particular forms of time dependence of the z/'s leading to desirable 
properties of the lognormal mixture dynamics and to a simple specification of the 
parameters in the model is carried out in Alexander and Brintalos (2003). 

A generalization of the mixture dynamics apparatus to multivariate underlying as- 
sets and possible applications to the pricing of basket options in presence of volatility 
smile are considered in Rapisarda (2002) and Brigo, Mercurio and Rapisarda (2002). 

3.4 Uncertain volatility geometric Brownian motion 

We conclude the paper by pointing out an important relationship between the lognormal- 
mixture diffusion dynamics and an analogous uncertain volatility model given by a 
geometric Brownian motion with uncertain volatility. 

In general it is known (Derman, Kani and Kamal (1997), Britten- Jones and Neu- 
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berger (2000), Gatheral (2001)) that every stochastic volatility model has a local 
volatility (i.e. scalar-diffusion) version that features the same marginal distributions 
in time (and thus the same initial prices for all plain vanilla options such as European 
calls). We may wonder whether our lognormal mixture diffusion dynamics ( TTOT) is the 
local volatility version of some stochastic volatility model. The answer is affirmative 
and we introduce the related model below. 

Consider the following uncertain volatility model: 

dS t = r(t)Stdt + t(t)S t dW t , S = s , (13) 

{(£ I— > V\{t)) with probability Ai, 
... 
[t i — ► f m (t)) with probability X m , 
£ independent of W and drawn at random at (an almost zero) time e > 0, 

where all z/'s are assumed to be regular enough and have a common value in [0,e], 
Vi{t) = v for all i and t < e. We assume £ to be independent of W and that the 
original probability space is large enough to allow for such a £ (otherwise we may 
define £ on a different space and then take the product space). Conditional on £, the 
process S is a geometric Brownian motion as in the Black-Scholes model. Thanks to 
independence, it is easy to show that for t > e > u, 

rn 

Q{S t eA\S u = y} = Y, \Q{St e A\S U = y,£ = u,} 
so that in particular 

Ps t \s u {x;y) = ^2kPs t \s u ,s(x;y,Vi), 

1=1 

i.e. the transition density of our uncertain volatility model is a mixture of lognormal 
transition densities, each corresponding to a volatility function v^. If we condition 
on an instant u > e, including the information on which value of £ has realized 
itself at time e < u, an information that is contained in the path of S up to u > 
e, the transition density between u and t reduces merely to a lognormal density 
characterized by the relevant z/j. 

By considering the case u — 0, it is immediate to see that this model has the 
same marginals as the lognormal mixture diffusion seen earlier, although it leads to 
an incomplete market. Then at the initial time it implies the same prices as the 
local volatility version seen earlier. Hedging is thus different and more complicated, 
and has to be based on additional hedging instruments. But what is interesting now 
is that transition densities are also known, not only marginals. This is not true for 
the local volatility version seen earlier. 

There is, however, a close relationship between the lognormal mixture diffusion 
dynamics ffTO]) and the uncertain volatility mixture dynamics ( TT3"j) . 
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Proposition 3.2. The lognormal mixture diffusion dynamics is the local volatil- 
ity version of the uncertain volatility mixture dynamics lp~3\) . The two models are 
linked by the relationship 

aUfx) = E{at) 2 \S t = x}. 
Proof. The proof is immediate by resorting to a variant of Bayes' formula: 



x 



E{m 2 \s t = x} = EiatfY, 1 ^ = Km = x\ = J2 E i^) 2 = ^}\ s < = 

fc=l fc=l 

m m 

= Y,Em 2 \St = x,£ = v k ]Q{Z = u k \S t = x} = J2"l(t)Qti = MSt = x} = <ZJt, 

k=l k=l 

since, by Bayes' formula, 

Q{£ = v k \S t = x} = Q{S t e dx\£ = v k }Q{£ = Vk}/Q{S t e dx} = A k (t, x). 

□ 

Remark 3.3. (Casting some light on the "zero terminal correlation" result of The- 
orem [37l\) The terminal correlation computed at time between the asset S and its 
average variance £ 2 (t)dt is easily seen to be zero, due to independence of £ and 
W . The same property is shared by the local volatility version, as was pointed out in 
Theorem \3.1[ Now we see that the local volatility version maintains the decorrelation 
pattern between volatility and underlying asset that is so natural for its uncertain 
volatility originator. The result of Theorem \3.1\ looks less surprising in the light of 
this result for the uncertain volatility version. 



3.5 Evolution of the volatility smile in the two models 

We now look at another important feature of the lognormal mixture diffusion that 
has not been investigated in our earlier papers. We are interested in looking at the 
smile evolution in time as implied by the model. 

We investigate this matter numerically as follows. We consider a set of parameters 
coming from a possible calibration of model (fTOj) with m = 2 to the foreign exchange 
market. The initial foreign exchange rate as from february 10, 2003 is So = 1.07 US 
Dollars for 1 Euro. Calibration of the model to market data provides us with the 
parameters (we assume the v's to be constant in time, except for a negligible initial 
time-interval [0,2e]) 

Ai = 0.9747, A 2 = 0.0253, v x (t) = 0.7572, u 2 (t) = 0.0899. 

The risk- neutral drift r(t) is taken consistently with the differences of interest rates 
in the domestic and foreign curves: indeed, we know that under the risk neutral 
measure, the drift of an exchange rate is the difference of instantaneous interest 
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rates between the domestic and the foreign markets, r(t) := Td{t) — rf(t), where 
interest rates are assumed to be deterministic and r^, 77, Rd, Rf denote respectively 
the domestic and foreign instantaneous interest rates and their integrals. The initial 
smile for an option maturing in one year, K > V(ly, K) is given in the Oy-vow of 
Table El We perform the following test. We consider options with one-year time to 
maturity T — t = ly set at future times t > 0, conditional on the average underlying 
being realized, i.e. conditional on S t = S t := E (S t ), the expectation taken under the 
risk neutral measure. We therefore price the option with our model ([TO]) , by resorting 
to an Euler scheme with time step of l/1000y (Antignani (2003)) and inverting then 
the corresponding Black Scholes formula. Indeed, let us define V(t, T, K) as the 
solution of the equation 

e- R ^E t [(S T - K) + \S t = S t ] = BSCall(4 K,T — t, R d (t, T), R f (t, T), V(t, T, K)). 

(14) 

By solving the above equation we infer what we call the conditional future smile at 
t for maturity T, i.e. K i— > V(t, T, K) / y/T — t. We will take t = ly,2y,3y,6y,7y 
and T = t + ly. We thus focus on options maturing in one year and on their smile 
as implied by model ffTOj) . which has been employed to compute the expectation on 
the left hand side of (fT4j) . Interest rates are taken from the market and assumed to 
be deterministic. Table [U reports the market values we used. We recall that in a 
deterministic interest-rates world R(t,T) = R(0,T) — R(0,t). 



T 


e -R d (o,T) 


e -R f (0,T) 


ly 


0.974454 


0.985738 


2y 


0.946724 


0.960891 


3y 


0.914757 


0.925555 


4y 


0.879548 


0.885228 


5y 


0.841922 


0.84227 


6y 


0.803363 


0.799019 


7y 


0.764915 


0.756466 



Table 1: Domestic and foreign discount factors (Euro and USD) 

We obtain the results reported in Table [2] as annualized percentage implied volatil- 
ities. We notice that the one-year smile flattens considerably in time (and similar 
considerations apply to longer maturity smiles). The initial smile (first row) shows 
an excursion from 16.17% to 10.72% and then up to 13.61%, whereas the last smile 
(7y/8y smile, last row) moves from 10.26% down to 9.55% and up again to 9.77%. 
As we see by looking at the rows of the table, the smile flattens considerably in time. 

Again, the one-dimensional diffusion model ([TO]) mimics the uncertain volatility 
model f[T3l of which it is the local volatility version. Indeed, since after time e we know 
the realization of £, we know which volatility V{ realized itself and, conditional on this, 
our process ([T3|) is a geometric Brownian motion and implies a flat smile. Therefore, 
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K/St 










t 


& 


0.8 


0.85 


0.9 


0.95 


1 


1.05 


1.1 


1.15 


1.2 


Oy 


1.07 


16.17 


13.73 


11.98 


11.02 


10.72 


10.84 


11.30 


12.26 


13.61 


ly 


1.08257 


11.23 


10.44 


10.10 


9.98 


9.96 


10.00 


10.10 


10.30 


10.62 


2y 


1.08628 


10.10 


9.95 


9.86 


9.83 


9.83 


9.88 


9.97 


10.07 


10.24 


3y 


1.08294 


10.60 


10.09 


9.92 


9.83 


9.81 


9.82 


9.85 


9.91 


10.02 


6y 


1.06401 


11.14 


10.20 


9.84 


9.71 


9.68 


9.66 


9.67 


9.72 


9.85 


7y 


1.05773 


10.26 


9.83 


9.67 


9.60 


9.57 


9.55 


9.57 


9.64 


9.77 



Table 2: Conditional future ly smile K t— > V(t, t + ly, K) for t = 0, 1, 2, 3, 6, 7 years 

after time e the smile flattens completely around the realized £ in the model ffTB"]) . 
and the one-dimensional diffusion version mimics this behaviour by progressively 
flattening the implied smile in time. 



4 Conclusions and further research 

In the present paper, we have found a candidate diffusion process whose marginal 
law follows a given evolving mixture of probability densities. We derived a stochastic 
differential equation (SDE) admitting a unique strong solution whose density evolves 
as a mixture of Gaussian densities. We introduced a seemingly paradoxical result on 
the comparison between the instantaneous and the terminal correlation between the 
obtained process and its time- averaged squared diffusion coefficient. As an applica- 
tion to option pricing, we considered diffusion processes whose marginal densities are 
mixtures of lognormal densities, showing how such processes can be used to model 
the market smile phenomenon. Furthermore, we have pointed out how the lognormal 
mixture dynamics is the one- dimensional diffusion version of a geometric Brown- 
ian motion with uncertain volatility, leading to an uncertain volatility Black-Scholes 
model, which allowed us to suitably reinterpret the earlier correlation result. Finally, 
we checked numerically the evolution of the smile in time and found that the diffusion 
model mimics again the uncertain volatility model by a substantial flattening of the 
smile implied curve in time. 
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